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Abstract. We study two exactly solvable five-dimensional thick brane world models in pure metric f(R) 
gravity. Working in the Einstein frame, we show that these solutions are stable against small linear pertur¬ 
bations, including the tensor, vector, and scalar modes. For both models, the corresponding gravitational 
zero mode is localized on the brane, which leads to the four-dimensional Newton’s law; while the massive 
modes are nonlocalized and only contribute a small correction to the Newton’s law at a large distance. 

PACS. ll.10.Lm Nonlinear or nonlocal theories and models - 11.27.+d Extended classical solutions; 
cosmic strings, domain walls, texture - 04.50.-h Higher-dimensional gravity and other theories of gravity 


1 Introduction 

The idea that our world might be a hyperspace (called 
brane world) embedded in higher-dimensional space-time 
(called the bulk) has been intensively considered in the 
passed two decades [TlEl[?lHll5l[51[71l51ini[IU] (for reviews, see 
IT2l ). This idea has changed our traditional knowledge to¬ 
ward extra dimensions. In early theories of extra dimen¬ 
sions, namely, the Kaluza-Klein type theories, the extra di¬ 
mensions are compacted to the Planck scale m While in 
brane world scenarios, depend on the model, the radiuses 
of extra dimensions can be as large as a few TeV' - 1 m 
ITU] , or several millimeters [5] , or even be infinitely large [5] 

0 . 

In one of Randall and Sundrum’s brane world sce¬ 
narios (the RS-II model) [5], the authors considered a 3- 
brane embedded in a five-dimensional anti de Sitter (AdS) 
space. Due to the nonfactorizable background geometry, 
the spectrum of four-dimensional gravitons of the RS-II 
model is consisted by a normalizable zero mode along with 
a continuum of non-localized massive KK modes. The nor¬ 
malizable zero mode corresponds to the four-dimensional 
massless graviton and leads to the Newton’s law. There is 
no mass gap between the zero mode and massive modes. 
So by intuition, the massive modes should cause a large 
correction to the four-dimensional Newton’s law. How¬ 
ever, after calculation Randall and Sundrum surprisingly 
found that the whole continuum of massive modes only 
contributes a small correction to the Newton’s law at a 
large distance [8]. In other words, the continuum modes 
are decoupled. 

In the set up of the RS-II model [8], the 3-brane has 
no thickness, and the geometry has a singularity at the 
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location of the brane. To evade this singularity, one can 
extend the RS-II model by replacing the original 3-brane 
by a smooth domain wall (called thick brane) generated 
by a background scalar field mm- Due to the config¬ 
uration of the domain wall, the bulk is not an AdSs space 
[TTJow. But the geometry is asymptotically AdS at the in¬ 
finity of the extra dimension. Thanks to this asymptotic 
behavior of the geometry, the gravitational zero mode is 
usually normalizable mmm- Besides, the authors of 
ref. [T6] found that at least for two mass points localized 
on the center of the thick brane, the continuum modes are 
decoupled provided the zero mode is normalizable. 

In addition to thick branes generated by scalar fields, 
there are also thick branes arise from pure geometry. For 
instance, by replacing the Riemannian geometry into a 
Weyl-integrable geometry, the authors of refs. mmm 
[2a constructed thick branes without introduce an addi¬ 
tional matter field. The normalization of the gravitational 
zero mode as well as the decoupling of the massive Kaluza- 
Klein (KK) modes are also studied therein. 

In this paper, we investigate another alternative for 
generating thick branes with only geometry. We assume 
that gravity is not described by general relativity, but by 
the so-called f(R) theories, where the Lagrangians are 
proportional to some functions of the scalar curvature 
R. The f(R) theories were created in the study of cos¬ 
mology mmm, and are mainly applied in cosmology 
nowadays (see and references therein). 

Nevertheless, there are works which devote to embedding 
branes, either thin or thick E3IM1EZ1 

into various types of f(R) gravities. 

Note that, all the thick /(i?)-branes considered in refs. [35l 
IMlIUniUHlIUUlHUlliTlI^lITU] are generated by a background 
scalar field. For f(R )~branes of this type, the tensor per- 
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turbation equation has been derived in ref. 0 a, but it is 
still unclear if these models are stable against the vector 
and especially, the scalar perturbations. To obtain reli¬ 
able thick /(R)-brane models, we must either proof that 
the solutions found in refs. 

are also stable against the vector and the scalar perturba¬ 
tions, or to find some new solutions whose stabilities are 
easier to proof. In this paper we adopt the second path. 

According to the well-known Barrow-Cotsakis theo¬ 
rem [23], a pure metric f{R) theory (referred to as the 
Jordan frame) is conformally equivalent to general relativ¬ 
ity minimally coupled with a single canonical scalar field 
(called the Einstein frame). This equivalence implies the 
possibility for constructing thick branes without introduc¬ 
ing additional matter fields. More importantly, it is much 
easier to analyze the linear stability of a solution in the 
Einstein frame. To the best of our knowledge, however, 
only refs. H5MS1 considered thick RS-II brane world solu¬ 
tions in pure f(R) gravity. In [55], the authors obtained 
a few numerical solutions. While the first analytical thick 
brane solution was reported recently in [46] . 

In this paper, we derive two analytical thick RS-II 
brane solutions in pure f(R) theories: one with a trian¬ 
gular f(R) and the other a polynomial f(R). The first 
solution is equivalent to the one of ref. @2)], despite an ap¬ 
parent difference. The second one is a new solution. These 
solutions will be presented in the next section. In section 
[3] we analyze the linear stability of these solutions in the 
Einstein frame by directly citing the results of refs. 03 
ITS] , Then in section Q] we show that the gravitational 
zero modes correspond to our solutions are normalizable, 
which implies that the four-dimensional Newton’s law can 
be reproduced on the branes. In section [5] by analyzing 
the asymptotic behavior of our solutions, we draw a con¬ 
clusion that for two mass points localized at the vicinity 
of the brane, the massive KK modes are decoupled, and 
only lead to small corrections the Newton’s law. In the 
last section, we summarize the main results of this paper. 


2 The model and the solution 

We consider pure f{R) gravity in five-dimensional space- 
time 


s = 2 ^| / d 5 xV^f(R), ( 1 ) 

where n 2 = 8ttG^ is the five-dimensional gravitational 
coupling constant, and g = det((/ mjv) is the determinant 
of the metric. In this paper we only consider the flat and 
static brane, for which the metric takes the following form: 

ds 2 = e 2A ^iy lv dx IJj dx 1 ' + dy 2 , (2) 

where e 2A ^ is the warp factor, ris the four-dimensional 
Minkowski metric, and y — x 4 denotes the extra dimen¬ 
sion. Throughout this paper, capital Latin letters M, N, ■ ■ ■ 
0,1,2, 3, 5 and Greek letters /u,, u, ■ ■ ■ = 0,1, 2, 3 are used 
to represent the bulk and brane indices, respectively. 


The Einstein equations read 

f{R) + Vr (4A 2 + i) - 6 f R A - 2 f R = 0, (3) 

and 

-8/ fl (i + i 2 ) + 8f R A-f(R) = 0, (4) 

where f R = df(R)/dR 1 and the over dots denote the 
derivatives with respect to y. By eliminating f(R), one 
immediately obtains the following equation: 

fR- fRA + 3f R A = 0. (5) 

For a specified A(y), eq. © is a second-order differen¬ 
tial equation for /r(j/). In case A(y) takes a simple math¬ 
ematical form, it is possible to solve f R (y) analytically. By 
inserting A and f R back into eq. Q, one can easily get 
the solution of f(R) as a function of y. Note that for the 
metric m, the scalar curvature R is related to y via the 
following equation 


R = —20A 2 - 8 A. (6) 

Once we get the expression of i?(y), it is not difficult to 
rewrite f R and f(R) as functions of R. 

Instead of starting with a simple /(i?), we prefer to 
begin with a simple A(y). For instance, we consider 

A = — nln(cosh(fcy)), (7) 

with n a dimensionless positive constant, and k another 
positive constant with the dimension of length inverse. It 
is convenient for us to introduce a dimensionless variable 
w = ky. In terms of w, the scalar curvature takes a simple 
form: 


R = 8nfc 2 sech 2 (iu) — 20n 2 k 2 tanh 2 (w), (8) 

from which we can express w in terms of R for an arbitrary 
n: 


w(R) = iarcsech 


v / 20n 2 + R/k 2 
y/8n + 20 n 2 


(9) 


The above equation makes it possible for us to get the 
analytical expression of f(R), at least for some special 
values of n. 


2.1 Case 1: n = 1, triangular f(R) 

We first consider the simplest case with n = 1, and 

A = — ln(cosh(w)). (10) 

Substituting eq. m into eq. ©, and only keep the 
=symmetric solution, one immediately obtains 


f R = cosh (a(w)) , 


( 11 ) 
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where the function a(w) is defined as 

a(w) = 2-\/3arctan ^tanh • (12) 

From eq. ©, one can easily obtain the solution of f(R(w)): 

f(R(w)) = 4/c 2 (3 — cosh(2w))sech 2 (w;) cosh (a(w)) 

— 8\/3fc 2 sech(t(>) tanh(tu) sinh (a(w )). (13) 


It is well-known that a pure f(R) gravity is confor¬ 
mally equivalent to a theory with a minimally coupled 
scalar in Einstein’s gravity [23] . The linear perturbations 
of the later case has been extensively investigated in lit¬ 
erature Thus, it is more convenient to discuss 

the stability of our solution in the Einstein frame. 

3.1 The Einstein frame 


Using eq. ©, we get 

f(R) = y (6 k 2 + 1?) cosh (a(w(R))) 

- ^ k2 ]J 4 80 - ^ ^ sinh (a(w(R))). (14) 

Note that f(R(w)) is an even function of w, so it makes no 
difference in choosing between the plus sign solution or the 
other in eq. ©. Therefore, eqs. ca and ra constitute 
the first analytically solvable f(R) brane model. 

Note that in ref. ESI. the authors also investigated 
thick RS-II brane solution in pure metric f(R) gravity 
with the same warp factor m- They also obtained an 
analytical expression of /(/?). Despite the difference in 
the mathematical expressions, it can be shown that both 
solutions are equivalent. The authors of [55] have shown 
that this solution is stable under tensor perturbation. In 
the next section, we will prove that this solution is also 
stable under scalar and vector perturbations. 


2.2 Case 2: n = 20, polynomial f(R) 

The second analytically solvable model appears when n = 
20 : 


A = — 201n(cosh(w)). (15) 

In this case, the symmetric solution of f R takes the form 

fn = 1 + 30 tanh 2 (w) + 65 tanh 4 (w ). (16) 

Using the same procedure as the last subsection, we obtain 
a simple polynomial solution 

f(R)=A + Cl R-^R 2 + ^R 3 , (17) 

where A = — 37 ygQ3 fc is the cosmological constant, while 

Cl = 2601 > C2 = — 41616 ’ c 3 = 39951360 are dimensionless 
constants. 

3 Linear perturbations and stability of the 
solutions 

In this section, we consider small metric perturbations 
around the solutions we obtained in the previous section. 
Our aim is to show that both of the solutions are stable 
against the metric perturbations to the linear order. 


Firstly, we define a new variable z, such that dz = e A dy. 
In terms of z, the metric can be written as 

gMN = e 2 A( ' z ' , r]MN- (18) 

Then, we introduce a conformal transformation 

gMN = ft{z) 2 gMN, (19) 

where I2(z) is a function of z. From now on, we will always 
use a tilde to denote a quantity in the Einstein frame. 
Obviously, gMN is conformally flat: 

gMN = d(z) 2 y M N- (20) 

Here a(z) = which will used in next subsection. 

Under the conformal transformation, the Ricci scalar 
transforms as m 

R = Q 2 R + 8g MN n{y M ^NR) - 20X, (21) 


where Vm is the covariant derivative defined by the con¬ 
formal metric gMN, and X = g AIN X n^- To con¬ 
tinue, let us first rewrite the original gravitational action 

CD as 


S = 


d 5 : 




f Ir „-5 


v2«i 


f2~ 5 R- Q~ 5 c 


( 22 ) 


where 


f R R-f(R) 

2«1 


(23) 


At this step, we only used the relation yf^g = Q~ b y/^-g. 
Next, we substitute eq. (1211) into eq. (1221) . and take f R = 
L? 3 , such that 


S = j d^x^iy^[R-l2n- 2 X}-^ 
By defining 


(24) 


■ = 2i/—y Ini?, V{(j)) = ft 5 d. 


(25) 


one can finally simplify the action as 

S = f *xV=§{^R-\x-V m . 


(26) 


This action describes a minimally coupled scalar field in 
Einstein’s gravity. The linearization of thick brane sys¬ 
tem with action (12H1) and metric (1201) has been thoroughly 
studied in [471149] , where the metric perturbations are clas¬ 
sified into tensor, vector, and scalar modes. Each type of 
these modes evolves independently, and none of the per¬ 
turbation equations relies on the explicit form of V ((/>). 
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3.2 Quadratical actions and stability 


Now we consider the linearization of a system defined by 
the action GUI) along with the metric GUI) . We need to con¬ 
sider perturbations comes from both the scalar field 4> and 
the metric < )mn , denoted by 8<f> and 8§mn = a 2 (z)h,MN, 
respectively. To obtain the equations for linear perturba¬ 
tions, one can expand the action GUI) to the second order 
of 8<f> and Jimn- The result can be found in refs. mm 
H8] , but here we use the one of ref. [ST] : 


S™ = 


+ 2«5 



d 5 xa 3 ^dMhNpd p h AIN — d M hd N H mn 
(d</>) 2 + 2h Mz (t>'d M 8(t) + <jJti8<j> 


- d M 8(j)d M 8(j) 


— ^ dphMN 


d p h MN 


i 


d p hdph 


+ 3 j{hd»h^ z -h zz h')}, (27) 

where d M = r] MN O n, h = r/ MN J imn, and the 

primes represent the derivatives with respect to z. Note 
that in this subsection, all the upper indices ji, v (or M, N) 
are raised by the Minkowski metric r] flv (or r] MN ). Note 
that for pure gravity around Minkowski background (a = 
1), S reduces to the well-known Fierz-Pauli action [52] . 

Following the procedures in ref. E5, we introduce the 
scalar-tensor-vector (STV) decompositions for the metric 
perturbatiorQ: 


K z = d li F + G ll , (28a) 

h^ u = T]^if + d^dyB + d^Cv + d v C^ + D^, (28b) 

where G^ are transverse vector perturbations: 

d ,l C,. = 0 = (29) 

and denotes the tensor perturbation, which is trans¬ 
verse and traceless (TT): 

ff'D lu , = 0 = D£. (30) 

The STV decomposition enables one to decompose S^ 
into three independent parts: 

SW = S& + S™ + S< 2 >. (31) 

Each type of perturbation evolves independently, and there¬ 
fore, can be analyzed separately. The vector and tensor 
sections are 

Sl 2) = \j d 5 xv»D^\, (32) 

5 t (2) = J + - ^i> ! '}(33) 

1 The STV decomposition method was firstly introduced 
in cosmology by Bardeen [53], and now is a widly accepted 
method in dealing with cosmological perturbations I54II55IIUU1 . 
This method can also be extended in the study of brane world 
perturbations I47II491I481I5T1 . 


correspondingly, where DU) = ri l ' lJ d IJ d l ,. The normal modes 
of the vector and the tensor perturbations are 

O'* = a*(G #t -Cp, = di, (34) 

respectively. 

The second-order action of scalar perturbations is more 
involved, it is composed by two parts I5H: si 2) = S s (2) + 

S s 2 2> where 

S; (2) = J d 5 xd 3 {^h zz - Zip' - 2k 2 <//<5</>}d( 4 V, (35) 

with %j} = F — i B', and 


c( 2 ) _ 
°s-2 — 


; J d 5 a;d 3 j — 3(/3D (4 V — 2>h zz \3^ip + Qip'ip' 


- 3-rh zz [h' zz + 4p') + 2re 2 8cj) + a' 


,d 2 V 

dej) 2 




2h zz (j)'8cj)' + <j>'{ti zs + 4ip')8<fi - (<5<(/) 2 ] }■ 


(36) 


( 2^ 

The variation SSg_{/8il> = 0 leads to the following con¬ 
straint equation 


3 jh zz - 3 p' - 2nlfiScp = 0. (37) 

Using this equation, one can eliminate h zz in the action 
m■ After a simplification, one finally obtains m 

sfj = J d 5 x g | □ + q" - ^ g J. (38) 

Here, g is a gauge invariant variable defined by 

g = ^a 3/2 (28 (j) (39) 

and 0 is a function defined as 

e = d 3 / 2 ^. (40) 

a' 

From quadratic actions m , ® and ([38]) . one can 
easily obtain the linear perturbation equations via the 
Hamiltonian variation principle 


8S^/5v^ = 0, (41) 

5S {2) /8D = 0, (42) 

ss^/se = 0, (43) 

and the hnal results are [5T] (see also S3): 

vector: □^ 4 ^u At = 0, (44) 

tensor: □ (4) l) /JJ/ + D " - ^ = 0, (45) 

a 2 

scalar: U^Q + g» _ 6 g = 0 . (46) 

U 
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Note that the tensor perturbation equation GSD has also 
been derived directly without using conformal transfor¬ 
mation by the present authors |44j . Obviously, the nor¬ 
mal mode of the vector perturbations has only zero mode. 
Therefore, our solutions are stable against the vector per¬ 
turbations. For the tensor and scalar modes, we introduce 
the following decompositions 


From eq. C%1) we know that the zero mode of the tensor 
perturbation takes the form 

Po{z) = A f T d 3/2 = A f T e* A j R , (54) 

where At is the normalization constant. The tensor zero 
mode po is normalizable provided 


Dp„(x x ,z) = e^e lpxxX p p (z), 

g{x x ,z) = e^ xX $ q (z), (47) 


where e pv is the TT polarization tensor. 

It is not difficult to show that p p (z) and $ q (z) satisfy 
the following equations 

AtA\pp = m 2 p p p , (48) 

A a A\$ q = M 2 <P q , (49) 

where m 2 p = -p p p p , M 2 = and 


At 


A s 


d . (giy 

dz + ~ a i ’ 
d 9' 
dz+ 9 


(50) 

(51) 


In the theory of supersymmetric quantum mechanics, the 
common structure of eqs. CHI) and CHI) ensures that both 
m 2 and M 2 are semi-positive definite, namely, m 2 , M 2 > 0 
for all p and q. Therefore, our solutions are also stable 
against the tensor and scalar perturbations. 

Note that in order to make the variables o, (f> and z 
well defined, our solutions must satisfy f R (y) > 0 for y £ 
(—oo, + 00 ). One can easily show that our solutions (flTT) 
and (USD satisfy this requirement. 


4 The normalization of the tensor zero mode 


Equation CHI) is in fact a Shrodinger-like equation 

- Pp + W(z)p p = rn 2 pp, (52) 

where the effective potential W(z) reads 

(a*)" _ (el A fl)" 


W{z) = 


~ 1 
a 2 

3aP_ 

4 a 2 


1 fs 
J E 


3 aT_ 
2 a 


3 o'f' R 


2 af R 


ul 

4 fl 


\fs (53) 


This expression consists with the result derived in ref. [44] . 

The spectrum of the tensor KK modes p p determines 
the effective four-dimensional gravity. Let us start with 
the zero mode po with Too = 0. A normalizable po leads 
to the four-dimensional Newton’s Law [81147] , Besides, the 
four-dimensional Planck constant is finite only when po is 
normalizable m- 


p+00 


1 = 


/ +00 

dzf R (z)e 3A ( z) 

-OO 


A/* 2 r+°° 

= ~jT J dwf R (w)e 


2 A(w) 


(55) 


Here we have used the relation dz = e~ A d^-. For both of 
our solutions, the above integration can be done analyti¬ 
cally. 

For the triangular model (fTTl) . the integration gives 


Ml 

2k 


cosh 



= 1 , 


(56) 


or 

M t ~ 1.953 Vk. (57) 

Similarly, for the polynomial model ED we get 


M t « 0.857^. 


(58) 


Thus, for both of our solutions, the gravitational zero 
mode is normalizable and can be localized on the brane, 
which results in the familiar Newton’s law on the brane. 


5 Correction to the Newton’s law 

To obtain an acceptable four-dimensional gravity, we have 
to require the massive modes p p with p > 0 do not lead 
to unacceptably large corrections to the four-dimensional 
Newton’s law (In this case, we also say that the massive 
modes are decoupled). For simplicity, we follow the study 
of ref. [TB] and only consider two massive points p\ and 
p ,2 located at z = 0. We denote the distance between pi 
and p, 2 as r. 

As have been addressed in ref. [16], the asymptotic 
behavior of the effective potential W(z) at \z\ 00 de¬ 

termines not only the localization of the zero mode, but 
also the decoupling of the massive modes. The localiza¬ 
tion of the zero mode requires that LF(|z| —> 00 ) > 0. 
If W(| 2 | —► 00 ) > 0, namely, there is a gap between the 
zero mode and the excited states, then we will obtain ex¬ 
ponentially suppressed corrections to the Newton’s law. 
The most interesting case is W(\z\ —> 00 ) —► 0. In this 
case, the scattering states start at m = 0, and the de¬ 
coupling of the massive modes becomes a delicate issue. 
An important result of ref. m states that if the potential 
W(z) —> P(j3 + 1 )/z 2 as \z\ —► 00 , the massive modes will 
contribute a correction AU oc 1/r 2 ^ to the Newton’s law 
at large distance (see also [57]). 

As depicted in figure [T] the effective potentials cor¬ 
responding to both of our solutions approach to zero as 
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W(z(y» 


y 



Fig. 1 . Plots of the effective potential W(z(y)). 

Wz 2 



Fig. 2. Plots of z 2 W(z) as a function of z. As \z\ — > oo, 
z 2 W(z) -> ¥. 


\y\ —> oo. They have the same asymptotic behavior in the 
z coordinate too, as the coordinate transformation is sim¬ 
ply a redefinition of y , and the shape of W will not change. 
Thus, our residual task is to prove that z 2 W(z) = /3(/3 + l) 
is a constant as |z| » 1, and to find out the exact values 
of P corresponding to our solutions. 

In fact, if z 2 W(z) approaches to a constant in the z co¬ 
ordinate, it should have the same asymptotic behavior in 
the w coordinate, in which all the quantities have analyt¬ 
ical forms. For instance, in the w coordinate, the effective 
potential is expressed as 


W(z(w)) = k 2 


s A ^)d w [e A ^d w {ei 



(59) 


and the variable z reads 


z(w) 



(60) 


For both of our solutions, z 2 W(z) can be analytically 
obtained. Instead of writing down the explicit expressions, 
we show in figure [2] that for both of our solutions 
lim^i^oo z 2 W(z) = ¥7 namely /? = 3/2. Thus, the cor¬ 
rections to the Newtonian potential AU oc 1/r 3 are sup¬ 
pressed at large r for both of our models. 


6 Conclusions 

In this paper, we studied two analytically solvable thick 
RS-II brane world models in pure metric f(R) gravity 
theories. Instead of starting with simple forms of f(R ), 
we began with two simple types of metric solutions, and 
derived the analytical forms of f{R). We obtained two 
types of f(R) gravities: a triangular one and a polynomial 
one. 

Then we studied the stability of our solutions against 
small linear perturbations, including tensor, vector, and 
scalar perturbations. We found that the solutions are sta¬ 
ble against all types of perturbations. 

In the end, we considered the reproduction of the four¬ 
dimensional Newtonian gravity. We first demonstrated that 
the tensor zero modes are normalizable and localized on 
the brane for both of our models; so the well-known New¬ 
ton’s gravitational law can be reproduced. Then we showed 
that for two static mass points localized at z = 0, all the 
massive tensor modes contribute a suppressed correction 
AU oc 1/r 3 to the Newton’s law. Therefore, the massive 
tensor modes are decoupled. 

This work compensates the studies of ref. gS] by adding 
a new analytical solution and offering a complete discus¬ 
sion on the stability of the solutions. Our procedures for 
finding analytical solutions might also be useful for cos- 
mologists. 
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